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AUSTRAL! 


This  report  sunn.irizcs  the  mean  wind  and  turbo 
lence  characteristics  in  and  above  a  junr.le-1  ike 
South  Carolina  coastal  forest.  The  data  reported 
are  from  one  trial  of  several  conducted  in  a  field 
stujy  by  Melpar,  Inc.  Very  little  interpretation 
of  the  data  is  given . 

The  vegetation  was  black  gun  -  red  maple.  data 
on  leaf  area  density  estimates  as  well  as  stem  densities 
are  tabulated.  The  mean  vector  wind  had  two  notable 
features:  a  nearly  constant  speed  in  the  lowest  two- 

thirds  of  the  canopy,  and  a  turn ing  of  the  wind  direction 
with  height  in  the  same  layer.  All  three  components  of 
turbulence  intensity  wore  larger  below  the  canopy, 
especially  t  lie  longitudinal  component.  Variance  spectra 
of  the  vector  wind  speed  indicated  that  an  inertial 
subrange  existed  at  all  heights  above  ant'  below  the 
canopy,  and  that  as  far  as  could  be  determined  the 
energy  dissipation  rate  '..’as  nearly  constant  with  height. 

The  diffusion  of  fluorescent  particles  from  a 
point  source  at  ground  level  snowed  that  the  initial 
drift  was  in  a  direction  60°  to  the  left  of  the  wind 
direction  above  tne  canopy,  but  beyond  a  distance  of 
ldb  n  from  the  source  the  centerline  of  the  plume 
gradually  turned  to  the  right.  The  lateral  cross- 
sections  of  jilurne  c  on  ce  n  t  r  a  t  ion  were  skewed  to  the  left. 
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ANALYSIS  OF  WIND  DATA  FROM  A  SOUTH  CAROLINA  COASTAL 

FOREST 

The  purpose  o  t  this  report  la  to  summarize  the  Been 
vine)  and  turbulence  chinctarlatles  In  and  above  a  jun  gl«- 
like  South  Caroline  coastal  forest.  The  data  reported  are 
from  a  field  study  conducted  by  Helper.  Inc.,  in  the  autur 
of  1964. 1  w*  have  performed  moat  of  the  wind  data  analyses 
ourselves  directly  from  raw  data  which  are  not  available  in 
published  fora.  It  is  intended  that  thla  report  should 
provide  selected  recent  observations  which  are  not  available 
elsewhere,  but  which  are  Important  to  the  advancement  of 
the  state  of  the  art  of  turbulent  diffusion  studies. 

We  will  select  one  particular  trial  (#11  ot  13  trials) 
during  which  the  equipment  was  most  reliable  and  the  speed* 
the  highest.  Wind  speed,  elevation,  and  azimuth  sensors 
(MRI  "vectorvanes")  were  used  at  eight  height  intervale  up 
to  60  meters.  The  forest  canopy  was  just  over  30  meters 
high.  The  vectorvanes  below  the  tree  crowns  were  operative 
although  not  very  active.  The  aystem  tolerances  for  the 
sensors  were  (1)  speed:  ±20  cm  sec'1,  (2)  azimuth:  ±8*. 

(3)  elevation:  ±3°,  The  data  ware  sampled  at  a  rate  of 
5.56  scans  per  second.  The  trial  period  vaa  135  minutes, 
beginning  at  2130  hours  EST,  i2  August  1964.  The  sit* 
was  under  strong  pre-»frontal  winds  from  a  cold  front  240 
km  to  the  north.  Upper  air  observation  ware  made  by  pilot 
balloon  and  theodolite.  The  mein  tower  site  was  300  metera 
from  a  forest  road  and  3  km  from  any  forest  border. 

Within  the  system  tolerances  for  temperature  Che 
temperature  profile,  in  the  layer  up  to  57  maters,  vaa 
for  all  purposes  isothermal  at  2 5  .  7 5  * C+0 . 0 5  * C  and  thus 
slightly  otable. 

1.  utScklP'fiuN  ui  Vegetation 


The  vegetation  within  a  150  meter  radius  of  the  main 
tower  site  was  relatively  homogeneous  in  composition  and 
stature.  Beyond  150  meters  there  was  to  the  vest  an  araa 
characterized  as  ’low  swamp",  and  to  the  east  an  area 
characterized  as  "pine  ridge".  These  areas  were  different 
in  composition  and  stature  from  the  tower  area.  It  can 
be  assumed,  however,  that  flow  transition  wss  not  seriously 
affected  by  this  difference  since  the  trea  tops  tended  to 


1.  A  summary  of  this  study  has  been  presented  by  Allison, 
Herrington,  and  Morton  (1). 


1 


■•if*  unlforaljr  In  the  horizontal.  The  tower  area  wee 
characterized  aa  "high  awaap".  Detailed  meaauremen ta  of 
the  foreet  physiognomy  and  composition  ware  aide  in  the 
high  and  lntaraadiate  awaap  araaa,  We  will  combine  theee 
data  for  purpoaee  of  characta r i zing  the  foreat  within  a 
100  meter  radiua  of  tha  tower  alte.  The  major  spaciea 
are  black  gum  and  rad  maple,  f-thsr  apeclea  are  tupelo 
gum,  bald  cyptese,  loblolly  pine,  poaiun  oak,  aweet  gum, 
and  holly.  Thaaa  "other"  treea  are  moatly  amall  under- 
atory  treea.  Tha  danaltlea  (>tcne  per  acre)  near  the 
tower  area  were  aa  followa,  indicating  the  difference  be¬ 
tween  the  high  and  the  intermediate  f.vamps. 


High  Swamp 

Intermediate  Swamp 

Black  Cum 

106 

145 

Rad  Haplo 

40 

130 

Other 

72 

33 

Total 

218 

309 

The  relationahip  between  diameter  at  brenat  height 
(dbh)  and  tree  height  v aa  not  a ign 1 f 1  can t ly  different 
for  the  high  and  intermediate  awaap  areaa.  On  the  baaia 
of  thia  aiailarity  in  atature  we  will  combine  the  two 
araaa  to  daacrlbe  tha  tower  area. 

Uaing  tha  following  aaauaptiona  and  atatlstical 
correlation* ,  it  la  poaaible  to  compute  a  vertical  dis¬ 
tribution  of  L(r)  ,  leaf  density  (surface  area  per  unit 
volume.)  from  detailed  phye  iog.ioml  c  measurements  of 
randomly  selected  red  maple  and  black  gum  trees,  it 
was  found  that  ragardlesa  of  tree  size  or  apeclea  the 
■•an  leaf  area  density  in  surface  area  (one  side)  per 
unit  tree  crown  volume  for  all  measured  crowns  was  0.41 
a- 1  with  an  error  of  0.07  m'l.  Since  thie  ia  a  relatively 
amall  error  considering  all  the  possible  factors  of  bio¬ 
logical  variability  we  shall  assume  that  this  density  in 
each  crown  was  a  constant,  and  that  the  effects  of  tree 
sice  and  species  ars  irrelevant  In  determining  how  leaves 
occupy  tha  available  apace  within  a  crown.  To  convert 


cross  sectional  area  at  a  given  height  level,  and  by  N, 
tha  number  of  trees  par  unit  ground  area  In  that  height 

class. 

Ihit  lit 

L ( z)“0 . 41N ( z) A( z)  m“l 
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The  density  of  trees  (sterna  per  ecre)wes  measured  for 
different  height  classes  and  we  use  this  for  N(z).  We  will 
have  to  make  use  of  none  statlBtlcal  correlations  to  obtain 
A(  z )  . 


The  average  cross  sectional  area  A(r)  can  be 
estimated  by  dividing  the  crown  voluses  D  of  a  given  height 
class  by  the  average  thickness  of  the  crowns.  D  was  not 
measured  on  a  height  class  basis.  The  non-linear  correla¬ 
tion  between  crown  volumes  D  ( to  3 )  and  trea  diameters,  dbh 
(Indies),  wee  as  follows! 

D-4. 5 (dbh) i*655 


This  correlation  was  fitted  to  the  combined  date  and  a 
standard  error  in  D  of  1972  was  obtained.  Although  the 
relative  error  is  large,  D  ranged  over  nearly  two  orders 
of  magnitude.  Considering  the  problems  in  sampling  and 
measuring  crown  volumes  we  shall  aasuma  that  all  the 
errors  were  sampling  errors,  and  that  the  above  correlation 
la  the  best  possible  estimate. 

Let  us  now  show  that  the  dbh  size  classes  can  be  con¬ 
verted  to  height  size  classes  with  the  aid  of  nonlinear 
statistical  correlation  formula  betvasn  dbh  (laches)  and 
tree  height  H(feet)  as  follows! 

H-13. l(dbh)°‘ 7 


This  correlation  was  fitted  to  the  combined  high  and 
intermediate  swamp  data  for  a  standard  error  In  H  of  29X. 

Uy  applying  the  above  assumptions  and  correlations 
to  the  tree  measurements  w«  were  able  to  obtain  a  distri¬ 
bution  of  leaf  area  density  L(z)  as  well  as  the  height 
Integrated  index  L^(z)  defined  as  follows: 

Lj(a)  “  J  Ldz  (d imenelon  lees ) 

It  should  be  noted  that  the  value  of  Lj  obtained  at  the 
top  of  the  forest  la  equivalent  to  the  "leaf  area  indax" 
(LAI)  used  In  agricultural  work.  LAI  la  defined  es  the 
total  leaf  area  of  the  crop  divided  by  Che  lend  area 
occupied  by  the  crop.  The  estimated  leaf  area  summary  Is 
given  in  Table  1,  and  is  shown  graphically  in  the  right 
sides  of  Figures  1,  3,  and  4. 

The  leaf  area  density  distribution  for  the  coastal 
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Table  1.  Leaf  area  distribution.  South  Carolina 
black  gum  -  red  maple  coastal  forest 


Height 

Leaf  Area  Density,  L(z) 

Height  Integrated  Index,  Lj 

m 

(m2/m3) 

( dimensionless ) 

3-55 

0.016 

0.028 

6,60 

0.047 

0.126 

9.65 

0.068 

0.303 

12.20 

0.026 

0.470 

15.lv 

O.O36 

O.56I 

16.35 

0.035 

0.603 

19- 4<i 

0.066 

0.760 

22. 45 

0.11U 

I.034 

25-50 

0.160 

1.451 

28.55 

0.045 

^.762 

31.60 

0.008 

1.044 

32.00 

0.000 

1.846 
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Table  2.  Tree  densities  versus  dimrteter  size  classes,  South 
Carolina  black  gum  -  red  majde  coastal  forest. 

diameter  at  breast  height  (dbh),  tree  densities 

size  class  (inches)  (stems  per  acre) 


3 

73- ^ 

4 

35.7 

5 

18.3 

6 

19-5 

8 

22.5 

10 

14.6 

12 

19.6 

14 

16.8 

16 

16.8 

18 

16.5 

20 

10.1 

i 


i 

? 


l 
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forest  had  two  local  maxima,  one  near  10  Mteri,  presum- 
ab ly  an  undsrstcry,  and  ona  near  25  meters,  the  auparlot 
crown  danaity.  Tha  LAI  for  the  foreat  wai  1.85  or  approxi¬ 
mately  oua-half  what  ona  would  expect  for  a  denaa  corn 
crop.  It  ahould  be  mentioned  that  etem  area  densities 
ware  determined  to  be  5.5Z  of  tha  laaf  area  densities. 

Other  data  which  are  useful  in  describing  the  vege¬ 
tation  are  the  summarized  tree  densities  (stems  per  acre) 
versus  dbh  (inches) .  These  data  show  that  in  the  tower 
area  tha  trees  are  "all-aged"  with  nearly  uniform  densities 
between  tha  sixa  claases  of  5  to  18  inches  dbh.  The  ab¬ 
sence  of  significant  variations  in  the  densities  in  this 
size  range  indicates  that  the  forest  has  been  undisturbed 
during  its  history.  See  Table  2. 

2.  MEAN  VECTOR  WIHD 

A  long-term  average  was  computed  for  che  vectorvane 
data.  Over  the  135-mlnute  sampling  period  these  computa¬ 
tions  usually  included  a  total  of  more  than  45000  data 
points  for  each  height  level.  The  azimuth  angles  were 
averaged  using  a  special  process  to  effectively  remove 
the  discontinuity  between  360*  end  0*.  The  elevation 
angles  averaged  to  zero  within  the  system  tolerances 
except  for  the  vectorvane  at  the  next  lowest  height  (9.2m) 
which  averaged  to  -13.4*  for  an  average  vertical  speed 
of  -2.3  cm/sec.  This  is  likely  to  be  an  instrument  mal¬ 
function,  however. 

It  ahould  be  noted  that  the  vectorvane  at  28.7m  was 
on  a  different  tower  and  in  a  strict  sense  lta  data  should 
be  excluded  from  the  analyses.  The  average  azimuth  angles 
at  59  m  and  at  3.8  m  could  possibly  be  erroneous.  We  be¬ 
lieve  tha  azimuth  angles  are  erroneous  at  59  m  and  3.8  m 
but  admit  to  a  preconceived  notion  of  what  should  occur. 

The  mean  azimuth  angle  at  59  m  is  smaller  than  values 
below  it  and  above  it  in  a  layer  where  this  should  not 
occur.  This  could  be  due  to  improper  referencing.  On  the 
other  hand  the  mean  azimuth  angle  at  3.8  m  could  have  been 
biaaed  by  the  Inertia  of  the  instrument  and  periods  of 
inactivity  more  than  instruments  above.  The  FP  cloud 
drift  is  a  better  Indicator  of  low  level  wind  drift  than 
the  instrument  at  3-8  m.  According  to  Morton  (2), 
project  manager,  the  speed  measurement  at  9.2  m  is 
"suspect"  due  to  malfunction. 

The  spaed  and  direction  of  the  mean  vector  wind  at 
different  height  levels  are  given  in  Table  3-A.  Supple¬ 
mental  wind  speed  and  direction  data  are  given  in  Table 
3-B,  which  are  upper  sir  theodolite  (pilot  balloon)  ob- 
ssrvatione,  a  Bendix  "Aerovane"  at  58  m  (Insensitive  to 
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Table  3A.  Speed  and  direction  of  the  mean  vector  wind.  135 
minute  average,  "vectorvanes". 


\ 

• 

Height  (m) 

Speed 

(cm  8ec"’^')+20,0 

Azimuth  (°)+8.0 

59.0 

447 

228.2* 

{ 

44. 1 

396 

247.1 

JU 

i 

37.8 

339 

251.6 

l 

32.4 

259 

241.6 

%£ 

f 

28.7 

135** 

267.3** 

f 

f 

18.  8 

29 

234.7 

b 

£ 

£ 

9.2 

11* 

215.3 

\ 

i 

3.8 

29 

250.2* 

i 

s . 

Table  3B. 

Supplemental  wind  observations. 

\ 

Height  (m) 

Type 

Speed  (cm 

sec”"*")  Azimuth  (° 

*’ 

i 

1500 ,  *** 

Theodolite 

2560 

— 

I 

£ 

610. 

Theodol 1 te 

1850 

265.0 

*r 

£ 

% 

* 

58. 

Ae  rovane 

521 

239.0 

V 

1.5 

FP  cloud 

24 

215.0 

f 

£ 

i 

I 

*  Assumed 

erroneous , 

£ 

f 

i 

S 

**  Not  on  the 

***  Estimated, 

same  tower  as 

o  the  r 

data  . 

turbulence) ,  and  fluorescent  particle  ( FT )  cloud  trace  at 
1.5  b.  Figure  1  shows  graphically  that  the  speed  of  the 
aaan  vector  wind  dropped  off  rapidly  through  the  tree  crowns 
and  was  nearly  constant  with  height  beneath  the  crowns. 

Also  Included  on  Figure  1  is  the  leaf  density  distribution 
L(x>  which  displays  the  principal  drag-producing  character¬ 
istic  of  the  forest.  Figure  2  depicts  a  turning  of  the 
wind  30*  to  the  right  froei  ground  level  up  to  the  tree 
tops,  which  is  verified  by  the  FP  trace  but  ve  have  ignor¬ 
ed  the  3.8  m  azimuth  measurement.  This  phenomenon  of 
turning  of  the  wind  has  been  verified  in  other  similar 
studies  not  reported  here  (1). 

Some  other  parameters  of  Importance  which  we  report 
here  ere  related  to  the  mean  vector  wind  above  the  tree 
crowns.  By  least-squares  fitting  of  the  logarithmic 
velocity  profile  to  the  speed  values  above  the  tree  tops 
we  obtained  the  following  mi crooe tecrologlcal  parameters 
and  their  vtandari  errors: 

1.  zaro  place  displacement,  2908  >80  cm 

2.  roughnesu  length,  z0,  16.81  ±8.90  cm 

3.  friction  velocity,  U^t  34.67±3.58  cm/sec 

The  fitted  values  of  speed  for  the  log  profile  were 
as  follows: 


HEIGHT  (m) 

59.0 

44.1 

37.8 

32.4 

3.  TURBULENCE  CHARACTERISTICS 


CALCULATED  SPEED  (cm/sec) 

449.1 

389.3 

342.2 

258.4 


The  high  rate  of  sampling  enabled  computation  in 
three  dimensions  of  the  intensity  of  turbulence  and  the 
coefficients  of  anisotropy.  Let  us  define  a  turbulent 
velocity  fluctuation  which  has  a  longitudinal  component 
u'  in  the  direction  of  the  mean  vector  wind,  a  lateral 
component  v'  that  lies  in  rhi  horizontal  plane  with  a  dir¬ 
ection  perpendicular  to  the  mea..  vector  wind,  and  a  ver¬ 
tical  component  w1.  The  sta..da.*d  deviation,  S,  of  each 
velocity  component  is  defined  as  follows: 
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Figure  2.  The  horizontal  direction  of  the  mean  vector 
wind  showing  turning  with  height  under  the  canopy.  FP  cloud 
drift  is  indicated  by  (x). 
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where  Che  bar  denote#  a  long-tern  averaging  process.  The 
fluctuations  are  lne tantaneoj#  variation#  from  the  long¬ 
term  components  of  the  mean  vector  wind. 

The  intensity  of  turbulence,  I,  in  each  component  direction 
is  defined  an  follow#! 


Iu“Su/v*  Iv-Sv/V,  Iw-Sw/V 

where  V  is  the  speed  of  the  mean  vector  wind. 

The  coefficients  of  anisotropy  Indicate  the  degree 
of  fluctuation  in  any  component  direction  relative  to  the 
fluctuation  in  the  longitudinal  direction.  In  the  lateral 
direction  the  coefficient  of  anisotropy  1#  Sv/Su,  and  in 
the  vertical,  Sw/Su.  These  express  the  relative  ehape  of 
the  turbulent  eddies. 

We  found  that  there  was  a  slight  temporal  trend  in  the 
wind  direction  azimuth  during  the  135-minute  averaging 
time.  For  that  reason  Sv  was  computed  by  averaging  27 
five-minute  values  of  Sv.  This  technique,  however,  filters 
out  all  variations  of  periodicity  greater  than  five  min¬ 
utes  at  the  same  time  it  eliminates  trends.  Alternative 
methods  of  removing  the  trends  were  tried  but  rejected  on 
the  basis  of  the  degree  of  subjectivity  in  a  curve-fitting 
process.  No  trends  were  removed  in  either  the  Su  or  Sw 
computation. 

The  longitudinal  component  of  turbulence  Intensity 
attained  a  value  greater  than  1.0  at  the  lowest  height 
levels  of  the  forest.  This  means  the  fluctuations  were 
usually  greater  than  the  mean  speed.  A  secondary  maximum 
occurred  just  above  tree-top  height  undoubtedly  due  to 
the  local  influence  of  tree  spires  on  the  flow,  see 
Figure  3. 

The  lateral  component  of  turbulence  intensity  also 
had  a  maximum  in  the  region  of  the  crowns.  The  vertical 
component  of  turbulence  intensity  was  distributed  parallel 
to  the  lateral  component,  except  the  values  were  cna- 
half  as  large  as  theme  of  t’r.a  later*.!  component.  All  Che 
data  of  turbulence  intensity  are  presented  in  Table  4  and 
are  shown  graphically  in  Figure  , 

The  coefficients  of  anlectopy  (Figure  4}  show  that 
above  the  tree  crowns  the  eddy  dimension  a  are  about  the 
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Table  4.  Turbulence  Intensity  components 


height 

(®) 

vv 

S  /V 

V 

VV 

59.0 

0.281 

0.330 

0.163 

44.1 

0.305 

0.335 

0.174 

37.8 

0.364 

0.347 

0.182 

32.4 

0.699 

0.387 

0.287 

18.8 

0.895 

0.445 

0. 32  7 

9.2 

1.222 

T  .  329 

0.191 

3.8 

0.  76  7 

0.329 

0.064 
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Figure  3.  Profile  of  the  three  u,v,w  components  of 
turbulence  intensity.  Shown  also  is  leaf  are  a  density 
L. 


CwErriCItNT  OF  Afijay  i  huKt  LEAF  AREA  DENSITY 

(nf1) 


Figure  4.  Profiles  of  tUe  coefficients  of  anisotropy. 
Also  shown  is  leaf  area  density  L. 
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lane  (isotropic)  for  both  components  In  the  horizontel 
plane  x,y.  But  the  vertical  component  shows  2 1 1  compress¬ 
ion  in  the  vertical  plane  (disc-shaped).  On  the  other 
hand,  the  eddy  shape  below  the  tree  crowns  has  a  dominant 
longitudinal  axis  but  is  compressed  about  2il  in  the  lat¬ 
eral  and  vertical  axes  (cigar-shaped).  These  data  are  not 
listed  since  they  can  be  computed  from  the  data  of  Table 
4 . 


It  should  be  emphasized  that  these  components  ere 
computed  in  a  locally  natural  co-ordinate  system.  That  is, 
x  is  in  che  direction  of  the  long-term  average  wind, 
which  is  not  the  same  at  each  height,  and  y  is  normal  to 
x  in  Che  horizontal  plane.  Thus  the  turbulence  intensities 
do  not  reflect  variances  due  to  systematic  turning  of  the 
wind  with  height. 

Because  of  the  high  sampling  rate,  the  turbulence 
data  lend  themselves  to  spectral  evaluation  of  the  var¬ 
iance  of  the  wind  speed  components.  Data  analysis  for 
the  variance  spectra  of  che  vector  wind  speed  follows, 

4.  VARIANCE  SPECTRA  OF  THE  VECTOR  WIND  SPEED. 


A  procedure  developed  by  Maynard,  McBride,  and  Stoner 
(3)  was  used  to  determine  the  variance  spectrum  of  the 
speed  of  the  vector  wind  (which  is  practically  equal  to 
the  speed  of  the  horizontal  wind).  This  procedure  defines 
the  variance  spectrum  as  the  Fourier  transform  of  a  "iven 
time  series  into  a  frequency  domain.  Following  Bend.it  and 
Piersol  (4),  the  variance  spectrum  was  written  an  the  real 
part  of  a  finite,  discrete  Fourier  transform,  and  a 
"Hanning  window"  was  used  as  a  smoothing  function  to  the 
spectral  estimates. 

The  number  f  spectral  estimates  is  determined  by 
the  number  of  lag  3  chosen  to  compute  the  autocorrelation 
function.  If  the  number  cf  lags  is  large  the  reliability 
of  each  estimate  Is  low  while  the  detail  of  the  spectral 
density  function  is  greater,  and  vice  versa.  We  chose 
120  lags  as  an  appropriate  number  to  obtain  sufficient 
detail  and  reliability.  The  raw  data  were  averaged  in 
non-overlapping  groups  of  5  before  spectral  processing 
(equivalent  to  a  0.90  sec  sample  in  real  time).  This 
reduced  the  number  of  data  points  to  9000. 

The  number  of  points  in  a  running  mean  was  determin¬ 
ed  by  tilal  and  error.  In  the  first  series  24,  121,  sc d 
201  point-  were  tried  and  it  was  found  that  24  points  gav* 
unbelievable,  unstable  results  while  121  and  201  points 
gave  Identical  results  at  high  frequencies.  At  low  fre¬ 
quencies,  however,  201  points  caused  an  apparent  loss  of 
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information.  In  the  second  series  the  number  of  points  In 
the  running  mean  was  varied  •;  o terns  tlcal  J.y  at  81,  101,  121, 
141,  161,  and  181  points,  and  tests  of  variance  were  made  as 
follows.  At  the  low  frequency  end  of  the  data  range,  it 
was  found  that  a  -5/3  power  law  was  approached  assympto- 
tically  in  the  spectral  density  estimate  for  all  height 
levels.  It  was  assumed  that  an  inertial  sub-range  was 
extant  at  the  low  frequency  end  of  the  data  range.  Thus 
the  spectral  dansity  function  was  expected  to  be  that  in¬ 
dicated  by  Lumlay  and  Panofaky  (6) 


Ek  -  0. 138  e2/3  k~5/3 


Where  e  ie  the  energy  dissipation  rate  and  k  is  wave  number 
(cycles  per  unit  length).  We  tested  the  effect  of  the  nua- 
ber  of  points  in  the  running  mean  by  finding  which  value 
gave  the  minimum  variance  of  e  in  the  low  frequency  (lew  k) 
end  of  the  data  range. 

Values  of  the  computed  energy  dissipation  rate  e, 
were  found  to  be  surprisingly  conservative  with  height 
through  the  forest  canopy  and  the  boundary  layer  above, 

The  computed  e  values  and  their  error  estimated  by  the 
minimixing  process  described  above  are  presented  in  Tsble 
5.  The  values  of  e  do  not  drop  off  rapidly  with  height 
in  the  boundary  layer  as  do  Ball's  values  in  Lumley  and 
Panofaky  (6).  However,  they  are  of  the  same  order  of 
magnitude  as  Bell  indicates  for  heights  of  10  to  100  m. 


The  feet  that 
through  the  forest 


e  is  practically  constant  with 
canopy  means  that  wherever  the 


heigh  t 
-5/3 
heigh  t 

in disting u is h&ule.  i u in  la 
plots  of  versus  k  where  considerable 


law  holds,  the  spectral  estimates 
level  a  will  be  practically 


at  different 


evident  in  the 

overlap  occurs  between  heights,  see  Figure  5.  L  It  should 
be  emphasized  that  these  data  were  taken  during  almost 
perfectly  isothermal  conditions  so  that  the  potential 


temperature 
stable) . 


gradient  was  positive  st  0.01  deg 


,-l 


(slightly 


At  higher  wave  numbers  (higher  frequencies)  the 
spectral  estimates  deviate  sharply  from  the  -5/3  power, 
especially  at  low  levele  within  the  forest  canopy.  These 


2.  For  convenience  of  plotting,  rpectral  estimates  at 
wave  numbers  greater  than  3x10”*  cycles  cm-^  were  averaged 
in  non-overlapping  groups  of  five. 
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Tabie  5.  Turbulent  energy  dissipation  rate,  e,  and  its  standard 
error,  Se ,  for  wave  numbers  less  than  the  upper  Unit  of  the 
inertial  subrange.  Below  the  canopy  (32*)  the  Halting  wave 
number  was  1.5xl0”3  cycles  ca“*. 


He igh  t 

e 

Se 

2  -3 

2 

a 

cm  aec 

ca  sec 

59.0 

158 

34.8 

44.  1 

120 

19.6 

37.8 

133 

44.4 

32  .4 

2  14 

83.0 

18.  8 

106 

41.2 

9.2 

149 

67.0 

3.  8 

54 

15.4 

Average 

134 

49.1 

i 

i 
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deviations  occur  over  a  significant  range  of  wave  nusbera 
that  they  cannot  be  assumed  to  be  an  artifact  of  data 
analysis.  In  fact  aliasing  of  the  data  is  not  suspected 
because  of  the  high  data  acquisition  rate  (5.56  scans  par 
sac)  and  because  the  original  data  were  averaged  in  non- 
overlapped  groups  <0.90  sac)  comparable  to  the  tiae  lag  of 
the  vectorvane  speed  sensor. 

the  deviation  from  a  -5/3  power  to  a  different  elope 
Indicates  the  limit  of  the  inertial  subrange.  At  the 
lover  height  levels  this  change  ia  due  to  momentum  ex¬ 
traction  by  the  drag  of  the  irregularly  distributed  rough¬ 
ness  elements  (leaves,  etc.).  The  critical  wave  number 
where  this  change  occurred  was  about  1s5x10“3  cycles  c»”^ 
for  the  18.8,  9.2,  and  3.8  a  height  levels,  ell  below 
the  canopy.  The  new  slope  of  the  energy  density  for 
wavs  numbers  greater  than  i.5xlO“3  cycles  cm~*  Indicates 
that  a  -3/10  power  of  k  defined  the  distribution,  see 
Figure  5. 

3.  DIFFUSION  CHARACTERISTICS 

Fluorescent  particles  (FF)  were  released  from  a 
bomblec  and  the  fallouc  was  measured  downstream  on  a  100  a 
horizontal  grid  network.  Different  colored  FP  were  used 
simultaneously.  The  FP  collectors  were  silicons  greased 
rods  callsd  "rotorods". 


A  brief  summary  of  the  diffusion  of  one  of  the  FP 
traces  released  at  the  height  of  1.5  m  will  be  presented. 
The  mean  size  of  the  FP  was  3.31  microns  with  a  size 
distribution  describsd  as  862  between  0.75  microns  and 
5.5  microns.  The  source  strength  wee  lO1^-  particles.  The 
settling  rat*  was  0*1  cm  sec”*.  The  total  count  of 
FP  at  the  1.5  height  was  determined  at  the  trial  conclu¬ 
sion.  The  centerline  oi  the  FP  plume  drifted  in  a  dir¬ 
ection  about  60s  to  the  left  of  the  wind  direction  above 
the  tree  crowns  for  s  distance  of  about  100  a  from  the 
source.  At  distances  greater  than  100  a  the  plume  center- 
line  gradually  turned  to  the  right  so  that  at  700  m  from 
the  source  the  direction  was  about  33*  to  the  left  of 
the  wind  direction  above  the  trees.  In  other  words  the 
plume  had  e  parabolic  trajectory  in  the  horizontal  x,y 
p  lane . 

The  distribution  of  FF  was  non-syaaet ric  about  the 
centerline.  Concentrations  were  skewed  to  the  left.  The 
maximum  total  count  (rotorods)  relative  to  the  source 
strength  decreased  along  the  centerline  distance,  x;  Che 
ratio  of  maximum  counts  to  the  source  strength  was  nearly 
proportional  to  l/x2.5#  Xhe  distribution  of  counts 
in  the  y-direction  perpendicular  to  the  centsrllne  is 
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shown  schematically  in  Figure  6  for  four  choean  di  *  ~ 
taacas  fro*  tha  source.  Tha  skewness  to  tha  laft  of  the 
pluaa  centerline  la  evident. 

6.  SUMMARY 

Vary  little  Intarpratation  of  tha  data  la  given  hare. 
It  vaa  fait  that  tha  praaantation  of  tha  data  la  tha  most 
important  iaaua  at  the  present  tine  in  order  that  investi¬ 
gators  nay  become  cogulxant  of  the  basic  wind  observations 
made  at  this  particular  site.  It  is  hoped  that  most  of 
tha  basic  questions  of  observation  have  bean  covered. 

Tha  vegetation  was  described  as  a  black  gum  -  red 
maple  coastal  forest.  The  leaf  area  density  distribution 
vaa  astimated  from  measurements  of  physiognomic  features 
of  the  vegetation.  We  indicated  that  the  forest  wae  un¬ 
disturbed  during  its  history  as  evident  in  the  even  dist¬ 
ribution  of  older  trees. 

The  mean  vector  wind  had  two  notable  features  in  this 
forest.  The  mean  speed  was  nearly  constant  up  to  heights 
about  two-thirds  the  height  of  the  canopy.  Secondly, 
the  mean  wind  direction  turned  with  height  about  30s  to 
the  right  in  the  same  layer. 

The  turbulence  intensity  changed  markedly  from  above 
to  below  the  forest  canopy.  All  three  intensity  components 
were  larger  below  the  canopy  and  the  longitudinal  compo¬ 
nent  was  especially  amplified. 

Variance  spectra  indicated  that  an  inertial  sub-range 
existed  at  all  height  levels  in  the  forest  for  wave  numbers 
less  than  1.5x10**^  cycles  cm~^«  There  was  no  significant 
difference  in  the  inertial  sub-range  energy  dissipation 
rate  above  or  below  the  canopy  within  the  errors  of  esti¬ 
mation  of  the  values.  At  the  lowest  heights  in  the  forest 
the  variance  spectra  had  a  -3/10  power  dependency  on  wave 
number  for  wave  numbers  greeter  than  1.5x10"^  cycles  cm~^. 

The  diffusion  of  fluorescent  particles  in  the  forest 
showed  that  initial  drift  of  a  source  at  ground  level  was 
in  a  direction  60*  to  the  left  of  the  wind  direction  above 
tha  canopy.  Beyond  a  distance  of  100  m  from  the  source 
the  centerline  of  the  plume  gradually  turned  to  Che  right 
and  approached  the  wind  direction  above  the  canopy.  The 
maximum  concentration  dropped  off  nearly  proportional  to 
l/x2*-’  and  the  lateral  distributions  were  skewed  to  Che 
left  of  the  centerline  due  to  turning  of  the  plume. 

More  analyses  will  be  performed  with  this  data.  A 
series  of  papers  are  in  preparation  interpreting  some  of 
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Figure  6,  Cross-sections  of  FP  concentrations  along 
the  plume  centerline  at  the  same  height  as  the  release 
height,  1.5a,  and  for  different  distances  from  the  source. 
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the  observations  in  Che  light  of  existing  theory 
and  In  comparison  with  the  results  of  other  workers. 
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This  report  summarizes  the  mean  wind  and  turbulence  characteristics  in  and 
above  a  jungle-like  South  Carolina  coastal  forest.  The  data  reported  are  from  one 
trial  of  several  in  a  field  study'  conducted  by  Melpar,  Inc.  Very  little  interpreta¬ 
tion  of  the  data  is  given. 

The  vegetation  was  black  guru-red  maple.  Data  on  leaf  area  density  estimates  as 
well  as  stem  densities  are  tabulated.  The  mean  vector  wind  had  two  notable  features! 
a  nearly  constant  speed  in  the  lowest  two-thirds  of  the  canopy,  and  a  turning  of  wind 
direction  with  height  in  the  sam;  layer.  All  three  components  of  the  turbulence  in¬ 
tensity  were  larger  below  the  canopy,  especially  the  longitudinal  component.  Variance 
spectra  of  the  vector  wind  speed  indicated  that  an  inertial  subrange  existed  at  all 
heights  above  and  below  the  canopy,  and  that  as  far  as  could  be  determined  the  energy 
dissipation  rate  was  nearly  constant  with  height. 

The  diffusion  of  fluorescent  particles  from  a  point  source  at  ground  level 
showed  that  the  initial  drift  was  in  a  direction  60°  to  the  left  of  the  wind  direction 
above  the  canopy,  bu’  beyond  a  distance  of  100m  from  the  source  the  centerline  of  the 
plane  gradually  turned  to  the  right.  The  lateral  cross-sections  of  plume  concentra¬ 
tion  were  skewed  to  the  left. 


UNCLASSIFIED 


